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四 元 数 Wishart 矩阵 分 布 及 其 在 通信 中 的 应 用 ” 


Æ dE dU 
(1- 烟台 大 学 数学 与 信息 科学 学 院 ， 烟 台 264005; 2- 华东 师范 大 学 数学 系 ， 上 海 200241) 
摘 要: 利用 奇异 的 四 元 数 Wishart 矩阵 的 密度 函数 以 及 特征 值 的 联合 密度 函数 计算 出 了 四 元 数 随 
机 MIMO 信道 模型 的 信道 容量 .此 模型 包括 信道 系数 矩阵 Hn, xn;， 它 可 分 为 nt S nw 和 ne > 
nr 这 两 种 情况 .文中 对 这 两 种 情形 分 别 计算 了 部 分 数据 . 
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1 引言 


关于 实数 随机 和 矩阵 的 多 元 统计 在 文献 [1] 中 已 有 详细 而 系统 的 研究 .在 对 复数 随机 矩阵 的 
WAF, Ratnarajah"!, Vaillancourt! 对 复 Wishart 随机 和 矩阵 (奇异 和 非 奇 异 )、 复 正 态 随机 
矩阵 的 密度 函数 和 复 Wishart 随机 和 矩阵 特征 值 的 联合 分 布 密度 函数 做 了 一 定 的 研究 ， 并 在 文 
BA [2] 中 给 出 了 复 随机 和 矩阵 在 信息 理论 中 的 重要 应 用 : MIMO Rayleigh 衰减 信道 的 容量 计算 . 
文献 [4-8] 在 四 元 数 随 机 和 矩阵 的 研究 基础 上 ， 建 立 了 四 元 数 的 随机 MIMO 信道 模型 ， 并 求 出 了 
该 信道 模型 的 容量 .以 下 是 要 用 到 的 文献 5,7] 中 关于 四 元 数 随机 和 矩阵 结论 . 

设 四 元 数 随机 矩阵 了 Y ~ QN. (H, In & X) (BY 服从 四 元 数 的 正 态 分 布 )， 其 密度 函数 为 


92mn 
q2mn | 5 | 2n 


TRABI W — Y"Y ~ QWa(n,Z)(n 2 m), BIW 满足 四 元 数 的 中 心 非 奇 异 Wishart 分 
布 ， 则 它 的 密度 函数 为 


exp (Retr( ~ 2771(Y — M)” (Y — M))). 


92mn 
QT, (2n)| n 


其 中 Retr(-) 表示 取 迹 的 实 部 ， 另 外 


exp (Retr(-2X1W))|W ?n7?m*1. 


QT, (o) = r”) I" Tla —2(j — 1), (o)»2(n- 1) 


j=l 
表示 四 元 数 和 矩阵 的 下 函数 . 文献 [7] 还 给 出 了 W 的 特征 值 的 联合 分 布 如 下 


92mn nem? —2m 


Qr. Ono, Gyr o P (C729 ^. Dp mn Ios — Aj) (4D), 


i<j 
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其 中 oFo(:) 是 四 元 数 超 几 何 函 数 ，D = diag(A ,Am)， 和 1,… ,Am 是 WW 的 特征 值 ， 当 2 = 
c?I WI, W 的 特征 值 的 联合 分 布 如 下 
92mn 42m? —2m 


QT (2m)QT m (2n)|a? 2n 





lppncsment Ilo: I3 AjJ* exp ( i, zo 5 X) {dD). (1) 
i—l 


i<j 
定理 1 EX Em x n 阶 的 四 元 数 第 阵 ， 秩 为 m(m 关中， 且 7(X) 是 关于 区 的 函数 ， 仅 通 
M XX" SX dx, M 
-T q mn aur 2 - 
jas dm Or, 2"! mM fW) W —WH o. 


证 明 HW = XXE =VAwVH (四 元 数 和 矩阵 的 谱 分 解 )， 其 中 Aw = diag (wi, , Wm) 
Wie ,wm 是 W 的 特征 值 ，V 是 广义 酉 矩 中 ， 则 由 文献 [8] 中 的 推论 1， 有 


dx 227" [|w (aw)(v7av) 


í-1 


及 文献 [6] 中 的 定理 6， 有 
9m42mn 


yHdy = 二 -一 一 一 ， 
(M QT, (2n) 


进而 ， 我 们 得 到 


/Co = fo f... 3X 


AW) WE mm" f yvHgav 
a Vini jr 


i 


gq?mn 
em Qr, eny LA 2m+1 f(W). 


2 FARNA Wishart E MIMO 中 的 应 用 


—^* MIMO 信道 系统 可 以 用 一 个 rr x n, 的 四 元 数 随机 矩阵 五 来 表示 ， 其 中 Ni, nr 分 别 是 
通信 系统 的 输入 、 输 出 数 . 该 通信 系统 如 图 1 所 示 ， 





图 1: 四 元 数 随 机 MIMGO 通信 系统 
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第 ;个 输出 接受 到 的 四 元 数 信号 可 以 写作 
yj = M higzi vj, 
i=1 


其 中 hj 是 第 i 个 输入 ， 第 ;个 输出 之 间 的 四 元 数 信道 系数 ，zx; 是 第 i 个 输入 的 四 元 数 信 
S. vj 是 服从 四 元 数 正 态 分 布 的 干扰 . 输出 端 接 收 到 的 信号 向 量 可 以 写作 


yi X1 vi 

hr >e Rna 
y2 . . . T2 v2 

一 . . . 十 

hin, Pes hn. 

Yn, Xn, Un, 
或 者 用 向 量 表示 如 下 
y= Hzrzcv, (2) 


AnyveQ"*, H e Qm, se Qxl, 输入 总 能 量 受 限于 p， 即 为 
E(zHz) < p， 或 trE(vzzP)«p. 


在 这 一 部 分 中 ， 我 们 将 通过 线性 模型 (2)， 求 得 四 元 数 MIMO 信道 模型 的 信道 容量 。 我 们 
假定 模型 是 块 衰落 模型 ， 且 在 多 个 独立 的 衰落 区 间 上 编码 ， 则 由 文献 [10]， 我 们 可 以 得 到 随 
机 MIMO 信道 的 Shannon 信道 容量 如 下 


C= Eg | log det (1s. t LH"H). 


其 中 万 ~ QN xn (0,1,, 8 X), W = HEH ~ QW, (nr, X). 上 式 中 的 对 数 函数 取 自 然 对 
数 时 ， 信 和 道 容量 是 以 nats 为 计量 单位 的 ， 车 对 数 函 数 取 以 2 为 底 的 ， 即 log。， 则 信道 容量 是 
以 bits 为 计量 单位 的 ， 本 文中 我 们 一 律 取 自然 对 数 . 

定理 2 ”如 果 输 入 能 量 限制 为 p, H QN. xn l0, In 9 X), Hn. 2 n:， 则 信道 容量 为 


ENT 2n,.—2n,4-1 p _1 
TEEF Spr log det (1,, + LW) exp (Retr(-25-1W))aW, 
其 中 W = HYH. 


证 明 通信 容量 为 
g= f logdet (1, + ^W), f(H)= te S AE: 
H Nt ni , q neni |x2n- 


由 定理 1 可知 
€ = me || logdet (Im + LW (Retr(-2x-! H7 H))dH 
= meum Jy og de (In. " ) exp etr(— )) 


" 22n«-n« p MNT 
i rm ja "EL (b; + 2w) exp (Retr(-2z- ! H7 H))dHdW 


22n«ni 


L 2n.—2n441 P | 9y—l 
GP. GS EP f, log det (In + £w) exp (Retr(-2x-! W))dw. 
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由 (1) 式 给 出 的 特征 值 的 联合 密度 ， 易 得 到 下 面 的 结果 . 
定理 3 ”如果 输入 能 量 限制 为 p, 五 ~ QNnxm(0,m 9X) Enr >ne Boe 


xf. oFy(-2X- p Aen Io. — A s (TEC db zX)) AP* 1 dg, 








izj 
其 中 . 
2n, n. sni 一 2nt 
K , 
QT, (2n;) QT, (2n; )| [277 
à >c > Àn, > 0 是 WW 的 特征 值 ， D = diag(Ai,:-- Ane]: 
我 们 称 忆 的 下 三 角 元 素 为 信道 的 相关 系数 ， 若 忆 的 下 三 角 元 素 是 非 零 的 ， 则 称 之 为 相关 信 
道 ， 反 之 ， 称 之 为 不 相关 信道 .下 面 我 们 将 举 一 个 数字 例子 ， 为 了 计算 方便 ， 我 们 取 n x 2 的 


不 相关 信道 ， 即 元 = o?J。. 我 们 可 以 得 到 nj. x 2 的 不 相关 信道 容量 如 下 


rm rug | eo (- ga) (oe EA) A 


(202) +2 (2n, 十 2) f? 
T(2n,)T (2n, — 2) J 


6(20?)?”r oo 


C -2 








exp ( 一 zz*) log (1 十 LA) sd 


exp (- zz ) log (1 十 LA) rd 


T(2n, = 2) Jo 
4(20?)?^-! T (2n, — 1) ?"--!T(2n, — 1) 2 " 
(OTOn)PQn.-2) Pu exp (7 grs) 10g (1+ EA) Ma 
| AQa?y?"-H T(2n, + 1) 5 1) 1 P "T 
"XGm)T(n.-—2) Jr (2n, 一 ] 9 (-z5*) log (1 i FAJA ^a]. 
o »" 上 述 数据 可 以 与 文献 11] 中 的 相关 数据 对 比 ， 我 们 可 以 发 


现 该 信道 模型 的 容量 远 远 大 于 同等 条 件 下 文献 [11] 中 的 复 Rayleigh 衰减 信道 容量 


























A1: ”四 元 数 随机 信道 容量 Hn, xni (nr > n) 

Nr 1dB 5dB 10 dB 15dB 20dB 

2 2.1380 x 10? 4.8343 x 10? 6.3114 x 10? 7.2384 x 10? 7.9182 x 10? 
4 2.2721 x 10? 4.5476 x 101? 5.6613 x 10? — 6.3310 x 10? — 6.8115 x 107? 
6 1.8790 x 101? 3.4835 x 107? 4.2342 x 10! — 4.6807 x 10!? 4.9995 x 10!? 
8 1.4231 x 107 2.5107 x 107" . 3.0085 x 107? 3.3030 x 107  . 3.5130 x 10?? 
10 1.0329 x10? 1.7590 x1025 2.0877 x 107 2.2815 x 107? 2.4194 x 10?5 
12  0.7314x 10? 12135 x10? 1.4295 x10?" 1.5586 x10? 1.6473 x 10?? 
14 0.5099 x 10? . 0.8286 x 1055 0.9705 x 105? . 1.0540 x10 — 1.1134 x 10?5 





3 “奇异 四 元 数 和 矩阵 在 MIMO 中 的 应 用 


在 这 一 部 分 中 ， 我 们 考虑 五 ~ QN. xn (0, In, 9 E), n, < nz 的 情况 ， 我 们 用 到 关于 奇异 
四 元 数 Wishart 矩阵 的 非 零 特征 值 联合 分 布 的 定理 自 ， 以 H ~ Q@N2xa(0, I5 8 X), X = o? H 
例 ， 部 分 计算 结果 如 表 2 所 示 (Wo? = 0.4). 
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表 2: ”四 元 数 随 机 信道 容量 Hn,xn (nr < m) 
Nr 1dB 5dB 10dB 15dB 20dB 
2  0.821x10!'  0.5282x10'? 0.6433 x 10? — 0.7117 x 10'? — 0.7606 x 10^ 
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Quaternion Wishart Matrix and its Applications in the 
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Abstract: The capacity of the quaternion MIMO channel model is computed by using the density 
function and the joint eigenvalue density function of the quaternion Wishart matrix. The model consists 
of the channel matrix 五 "xn with n; € n. and ne > nr. Some data are further computed for both 


kinds of matrices. 
Keywords: quaternion MIMO channel; channel capacity; distribution function of quaternion Wishart 


matrix 
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